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1 Introduction 

It is known that the response of a scintillation detector can't be approximated by a 
symmetric shape since the line skewness is not zero [l] (see also discussion below). 
An example of the situation where the deviations of the line shape from a gaussian 
can lead to systematic errors is the search for the effects on the tail of beta-spectra: 
smearing of the spectrum due to the detector's finite resolution provides a stronger 
underlying background in comparison to what one would expect in the case of a 
gaussian line shape. 

The purpose of this work is to provide a simple analytical expression for the 
asymmetrical shape approximating the corresponding ideal scintillation detector 
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response for average scintillation intensity counting from tens to hundreds of regis- 
tered photoelectrons. 

2 Ideal scintillation detector 

The statistical properties of a scintillation detector response were studied by Breit- 
enberg [T] and independently by Wright [2]. They showed that the relative variance 
Vq = of the scintillation detector pulse height is: 

VQ=VT + {l+VT){Vn--) + —^, (1) 

where vt is the relative variance of the photons transfer efficiency, /i is the mean 
signal registered at the photomultiplier (PMT) anode, measured in photoelectrons 
(p.e.), n is the mean number of photons produced in a scintillation event and Vn 
is a relative variance of the number of photons (which reduces to ^ in the case of 

the normal or Poisson variance), and Vi = (^^^ is a relative variance of the single 
photoelectron response (s.e.r.) of the photomultiplier (gi and ai are mean position 
and variance of the single p.e. peak). 

We will consider an ideal detector with the following features: 

1. anode signal for a single registered photoelectron is described by normal dis- 
tribution; 

2. the photoelectrons are registered statistically independent; 

3. the number of registered photoelectrons (p.e.) n for a monoenergetic source 
with a mean number of registered p.e. /i, follows a Poisson distribution, 

Pin) = ^e-^; 

4. intrinsic line-width of the scintillator is negligible, the variance of the number 
of scintillation photons is normal; 

5. the detector is spatially uniform, i.e. events with the same energy produce 
identical responses on the average at any position inside the detector; 

6. noises in the system are negligible. 



2 



As it will be shown below, condition (1) is essential only when registering on the 
average small numbers of p.e. in an event, /i < 8. Condition (2) is satisfied 
practically automatically in the case of detector with many PMTs working in single- 
electron regime, but could be questionable for a scintillator crystal coupled to a 
single PMT. Assumption (3) is natural, but (4) can need further validation in a 
real- world detector (especially in the case of the solid-state scintillators). Condition 
(5) is difficult to satisfy for large volume detectors, but in the case of a spatially 
non- uniform detector it is enough to introduce an additional parameter vt, defined 
above, to improve the fit quality. An example of fitting the ^''C beta- decay spectrum 
in a large volume non-uniform detector will be given below (see subsection I5.3p . 

In [3j the case of a real scintillation detector with many PMTs is considered, 
and it is shown that in the above assumptions ([T]) reduces to: 



where I)i is a relative variance of the single photoelectron response averaged 
over all PMTs of the detector. Thus the scintillation detector consisting of many 
identical PMTs, surrounding the scintillator can be considered as one PMT with an 
extended photocathode. For this reason the terms "PMT" and the "detector" will 
not be distinguished in the following discussion. 

If the PMT response (anode output pulse height q) to precisely n photoelectrons 
is fnil), and the number of the registered photoelectrons is distributed according 
to distribution P(n), then the PMT response function can be written as f{q) = 
'^P{n)fn- The PMT response function here is the probability density function 
(p.d.f.), it is normalized to the unity. At the absence of photoelectrons at the input 
of the electron multiplier (n = 0) the PMT is registering the noise of the system in 
accordance with the p.d.f. fo{q). Using the assumption of statistical independence 
of the registered photoelectrons one can write the p.d.f. of registering precisely 
n photoelectrons as a convolution of n independent single-photoelectron signals 
/n = /i ® ••• ® /i- If /i is described with a normal distribution, then follows a 
normal distribution as well, with mean n ■ qi and variance cr„ = y/nai. 

With a proper choice of fi{q) function the p.d.f. of the PMT response can be 
constructed at any mean scintillation intensity 

f{q) = Y,nn)Uq) = P{0)fo{q) + $^P(n)/„(g) ® /o(g). (3) 

n=0 n=l 
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The Fourier transform of (l3|) gives the characteristic function: 

Xis) = P(0)xo(s) + Yl Pin)xi{s)Xois), (4) 

n=l 

where Xii^) and Xoi^) are characteristic functions of the single photoelectron 
response and noise, respectively. 

For the case of the Poisson distribution of the probability to register precisely 
n p.e. in a scintillation event of mean intensity n p.e., the contributions from 
n = 1,2... p.e. can be summed in and ([4]) can be rewritten in a more compact way: 

X{s) = e-^xo{s) + J2 ^e-^xUs)Xo{s) = e^^^^^^^-'\o{s) . (5) 

n=l 

The analogous formula can be obtained for the generating function by using the 
elementary facts from the theory of branching processes [4]. In fact, omitting the 
noise term, equation ((51) corresponds to a 2-stage cascade device: the photocathode 
and electrostatic focusing system providing on the average fi Poisson-distributed 
photoelectrons at the entrance of the electron multiplier with generating function 
6*2 (s) = e'^*^'*"^-*; and the electron multiplier itself with a single photoelectron re- 
sponse at anode fi{q) with corresponding generating function Gi{s). The resulting 
generating function has the same form as (l5|): G{s) = G2{Gi{s)) = e^^'^'^^'^'>~^\ 
except of the noise term Xo{^)- 

Omitting the noise term, equation ((31) gets the form /(x) = X]n=o -^(^)/"(^) 
with characteristic function xi^) = e^^^'^^'^^~^\ which defines the so called com- 
pound Poisson distribution: the probability distribution of a "Poisson-distributed 
number" of independent identically-distributed random variables [5]. In our case 
the elementary distribution is the s.e.r., and the number of the independently reg- 
istered photoelectrons varies in accordance with Poisson distribution (assumptions 
2 and 3). 

The inverse transform of ([5]) in some special cases of Xi{^) can be performed 
analytically, for example, the case of an exponential single photoelectron response 
was considered by Prescott in [6]. 

An example of realistic function is shown in FiglH This is the average 

response observed for the ETL9351 photomultiplier used in the Borexino detector 
[7], the measured mean relative variance over a set of 2200 PMTs selected for the 
detector is Vi = 0.34 [8|. If the single photoelectron response of PMT and noise 
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function are known, then formula (JH]) can be used to construct the PMT response 
for any fi for which the basic assumptions are valid. The method based on the use 
of transform (l5l) has been successfully applied to fit the experimental spectra ob- 
tained with electrostatically focused hybrid photomultiplier tubes for few registered 
photoelectrons (/i = 2.66 and = 6.36 p.e.) in [9], where formula (l5|) was called 
"light spectra sum rule". 

It should be noted that single photoelectron spectra of the photomultiplier stud- 
ied in [9] has a very narrow single p.e. peak, so that the detector response to 
fi = 6.36 has "fine structure" peaks around the values corresponding to integer 
numbers of the registered charge. In this article we consider a case of fi^ fiQ with 
/io big enough to make the contribution of the first resolved n— fold photoelectron 
peaks to be negligibly small. The parameter /io can be obtained from the following 
considerations. The PMT response to precisely n p.e. (n-fold peak) with increase 
of n converges very fast to a normal distribution with q = nqi and = naf as it 
follows from the central limit theorem. In practice the PMT response to as low as 
n > 3 p.e. can be approximated by a gaussian, see i.e. [lO]. The (n-l)-fold and 
n-fold peaks are not resolved if the half width on the half heights resolution of the 
n-th peak is worse than ^qi. ^/2\n2 ^/naf > i.e. n > The contribution 
of responses from few photoelectrons decreases very fast with the increase of fi. It 
is easy to check that the condition P(0) + -P(l) + -P(2) < 0.01 is satisfied already 
at yUo — 8 P-G. In this case instead of the real shape fi{q) of the PMT single elec- 
tron response one can choose the gaussian approximation for the function 
with mean qi and variance ai coinciding with the corresponding parameters of the 
real-shape function. Indeed, the response functions for 3 and more p.e. are well 
approximated by a normal distribution, and 0,1 and 2 photoelectrons contribute 
less than 1% to the total spectrum (see also FiglSj). 

In such a way an ideal detector response is described by the inverse transform 
of (I5|) with Xi{^) corresponding to the characteristic function of a gaussian with the 
mean value and variance of the corresponding single photoelectron response: 

=e-5'^'^'e'^i^ (6) 

In the following discussion we call the "ideal" detector response obtained from 
([5]) by using Xi{s) from (I6|), and we let the "real" detector response to refer to 
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Figure 1: An example of the single electron response 



(I5j) with Xi{s) obtained by transforming the real shape of the single photoelectron 
response. The difference between the "real" and "ideal" scintillation response van- 
ishes very fast with the increase of /i (at /i > 8 p.e.). We have chosen the gaussian 
shape for s.e.r. for convenience, but any appropriate s.e.r. line shape can be used 
(with a relative variance that of real s.e.r.). This is illustrated in Figl2l where the 
theoretical photomultiplier responses for /i = 3 p.e. obtained for 3 different s.e.r. 
function (realistic from FigUl gaussian and rectangular) with the same mean value 
and variance, are plotted. One can see that the difference is noticeable only at the 
registered charge Q < 3 p.e., the tail of the PMT response is modeled equally good 
with the gaussian and rectangular s.e.r. functionJll- 

^So, attempts to evaluate the single electron response spectrum at /i > 1 seems to be senseless 
for the PMT spectra with unresolved s.e.r. {vi > 0.18), in the best case one can succeed to extract 
qi and wi values, but not the details of the s.e.r. shape. 
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Q, p.e. 



Figure 2: Photomultiplier response obtained for 3 different single electron response 
functions for the case fi = 3 p.e. 

3 The normal distribution limit case for ideal 
scintillation detector response 

The ideal detector response converges quickly to the normal distribution as /i grows. 
In fact, the Poisson distribution of the primary photoelectrons at the input of the 
electron multiplier converges to a normal distribution for big /i. The variance in 
the multiplication of the photoelectrons arriving at the electron multiplier, for high 
/i values can be considered roughly the same for all possible values of the registered 
number of photoelectrons (cr(/i+Ayu) = ^/7t(l+t'i) + |^^^AyU+.. ~ So, in the 

big /X limit the ideal response converges to the convolution of two gaussian processes 
which give a normal distribution with the mean value and variance, respectively: 

a' = il + v^)■^ = iql + al)^,, (7) 
/i 

coinciding with the values found above considering statistical properties of the scin- 
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tillation registration process. We assume that the scale is cahbrated in photoelec- 
trons, i.e. qi = I (otherwise it is necessary to pass to variable ^). The characteristic 
function for a gaussian p.d.f. is: 



1 „2 „2 



X{s) = e—^"^' e*"' (8) 

and it is apparently different from an ideal shape characteristic function ^ with 
Xi(s) from (jni). Moreover, one can calculate the moments of the ideal scintillator 
response from its generating function: 

Mn = (-)"^ 1.^0, (9) 

and check that only the first two moments of the gaussian and ideal responses are 
equal. The third central moment calculated for the ideal response is M| = (l+Sf 
which neither coincides with that of a normal distribution (it is simply zero), nor 
converges to it with increasing /i. Only the skew s = — which is a measure of the 
distribution asymmetry, indeed converges to zero slowly enough as 

Although the normal approximation of the scintillation line shape is quite com- 
mon [T], there are situations in which its use leads to systematic errors in the 
parameter definition. Two examples will be considered below (see section [5|). In 
order to resolve this problem, a better approximation of an ideal scintillation shape 
is needed. 



4 The generalized gamma distribution as a limiting 
case for the ideal response 

We will search for a function with the following properties: 

1. the function converges to a normal distribution for — »• oo; 

2. it has the mean value and variance coinciding with that of the ideal scintillator 
response; 

3. it approximates the ideal scintillator response better than a conventional nor- 
mal distribution; 

4. it is asymmetric with a skew decreasing as and gives a better approxima- 
tion of the distribution tail. 
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In literature the successful usage of the 2-parameter gamma- distribution to ap- 
proximate the output pulse height spectra of scintillation detectors is reported, 
with better results in comparison with a normal approximation [II],[l2]- We were 
not able to get a good agreement with the response function of an ideal detector 
using the above- mentioned distribution, so we have chosen a power transformed 
gamma distribution (also known as generalized gamma distribution) as a candidate: 

f{x;m,a,f3) = -^/3™"a;™"-^e-(^^)'". (10) 

The distribution describes a variety of well-known 1 and 2-parameter probability 
laws as special cases; more details regarding the distribution properties can be found 
in [13]. A physical basis for the generalized gamma distribution has been discussed 
by Lienhard and Meyer in |14| . 

We start by fitting the ideal scintillator response for different /i values using 
([Toll with 3 free parameters. It has been discovered that over a wide region of fj, the 
value of parameter m is close to 2, thus we fix it at this value and use the following 
distribution as an approximation of the ideal shape response (redefining from 

i^u^ as py. 

g{q;a,(3) = 2(3^r~\a)q"^-'e-^''\ (11) 

with parameters a and /3 providing equality of the mean value and variance of 
( ITTll to the corresponding values of the ideal scintillation response. It is easy to 
check that the moment of order n of the distribution (fTTll is: 



r(a + f) 



T{a) • 

The parameters a and P can be defined from the system of equations: 

^-f'- na) P ^2) 

A recipe for the approximate solution of the system is given in Appendix A. An 
alternative way of calculating the parameters a and [3 based on the equality of the 
first two even moments of (fTB to the corresponding values of the ideal scintillation 
response, is presented in Appendix B. 
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It is important to stress that a special case m = 2 is found in many physical 
applications: in hydrology it is known either as hydrograh distribution [14], or in 
countries where the Russian hydrology school has become more familiar, as the 
Kritskiy- Menkel distribution |15| : in radio-engineering variants of the generalized 
gamma-distribution are widely used to describe radio waves propagation in fading 
environment (Nakagami distribution [IB]); some further examples can be found in 

[ni- 

In the limit a ^ oo the distribution g{q) converges to a normal distribution 
[T8] . the condition 2 is satisfied automatically, conditions 3 and 4 have been checked 
numerically in a wide range of /i values. As it can be seen in FigjS] the generalized 
gamma distribution approximates the ideal response better than a gaussian. FigH] 
presents results of numerical calculations of the deviation of the gaussian (with the 
mean value and variance that of an ideal response) and the shape obtained with 
([TTll from the ideal response calculated as: 

/ \9{q) - f{q)\dq, (13) 

and has a simple mathematical interpretation. In FigH] one can see that the devi- 
ation of the generalized gamma-distribution from the ideal one calculated by using 
(fTS]) is an order of magnitude lower than that in the gaussian distribution case. 

The quality of the fit in the tail has been checked by calculating the integral 
in the region [/i + 2a; oo] for the ideal and generalized gamma- distributions. The 
integral of the gaussian in this region is constant defined by the complementary 
errors function: 0.5erfc(-v/2). The cumulative distribution corresponding to the 
density f fTTl) is: 

G{x) = / g{x)dx = 'y{a, f3x^), (14) 
Jo 

where 7(0?, x) is the normalized incomplete gamma function. The integral in the 
tail is 1 - G(/i + 2cr). 

Integral in the tail for the ideal response was calculated by using the original 
definition ([3]): 

^In [13] the case m = 2 is called Stratonovich distribution. We were unable to find the 
corresponding reference in literature. 
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Figure 3: Comparison of the ideal scintillation response with the gaussian and the 
model by means of a generalized gamma distribution for /i =10, 20, 50 and 200 
p.e. Responses obtained by using the realistic s.e.r. function (see FigH]) are not 
distinguishable from the ideal scintillation response in all the above plots. 
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Figure 4: The deviation of the response constructed by using the generalized gamma 
function from the ideal one is an order of magnitude lower than that for the corre- 
sponding gaussian. The deviation was calculated by means of (fT3l l. 



with Nmin = max{[fi — 2cr],0) and Nmax = + 5cr. The results are presented in 
FiglHl One can see that the gamma distribution gives a better approximation of 
the distribution tail than the gaussian one. 



The most probable value of distribution ( ITTi l corresponds to q = y ^{<^ — ^) 

[TT] . it can be seen that q is shifted to the left from the mean value by ~ ^-j^- 

5 Two examples 

The precision of the description of the spectra of a real scintillation detector with 
respect to different approximations of the response function has been verified by 
using both the real data of the Counting Test Facility (CTF, [H]) of the Borexino 
detector [7], and the data obtained with the Monte Carlo model of the CTF detector. 




n=N, 
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Figure 5: The ideal scintillation response tail is reproduced very well for yU > 8 p.e. 
The corresponding gaussian response tail does not depend on /i and is defined by 
ierfc(v^). 

In the present article we consider only the MC data, the results of comparison of 
the theoretical model with the real CTF data will be presented by the Borexino 
collaboration. 

The large volume liquid scintillator detector CTF is a prototype of the solar 
neutrino detector Borexino. The CTF was used to develop the methods of deep pu- 
rification of the liquid scintillator and water from the natural radioactive impurities. 
The CTF consists of 3.7 tones of liquid scintillator on the base of pseudocumene 
(C9H12), contained in a transparent spherical inner vessel with a radius of 1 m, and 
viewed by 100 photomultipliers (PMTs) mounted on an open spherical steel support 
structure. The PMTs are equipped with light concentrator cones to increase the 
light collection efficiency; the total geometrical coverage of the system is 21%. The 
radius of the sphere passing through the opening of the light cones is 2.73 m. The 
entire detector is placed inside a cylindrical tank with water, which provides shield- 
ing against external gammas. On the bottom of the tank another 16 PMTs are 
mounted to identify cosmic muons by their Cherenkov light produced in the water. 
The detailed description of the CTF detector can be found in [19]. The CTF has 
been in operation since 1993. At present it is in its third data-taking campaign 
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(CTF3) with the main goal of tuning the purification strategy for the Borexino 
detector. The data collected with an upgraded version of the CTF were used by 
Borexino collaboration in order to search for a number of possible manifestations 
of non-standard physics, a review of experimental results can be found in [20] . 

The Monte Carlo model of the CTF detector was developed on the basis of 
EGS-4 code [2T] to check the validity of the background interpretation. It accounts 
for the dependence of the light yield on the energy (ionization quenching) and on 
the position where energy was deposited inside the detector. The model has been 
calibrated with the CTF data and describes the CTF experimental spectra with 
a satisfactory precision. For the purposes of the present work, the model of the 
detector response was changed to take into account the deviations of the response 
function from the normal one (the standard program uses the normal approximation 
of the response function). 

5.1 Monoenergetic line 

The detector response to the monoenergetic particle has been modeled with the 
MC method. The particle energy was chosen in order to provide the number of 
registered photoelectrons, fi = 150 p.e. The number is big enough to ensure good 
approximation with a gaussian shape. Indeed, the processing of the CTF data by 
using this approximation was successfully applied even for lower values of the mean 
registered charge [22] . 

The response of the detector was generated in the following way. First, the mean 
number fj,Q of p.e. registered at one PMT was defined as /xq = where Npmt 

is the total number of the PMTs in the detector. Then in each event for each PMT 
the Poisson- distributed number K of registered p.e. was generated, and, finally, 
the registered anode charge was simulated using the gaussian approximation of the 
PMT signal with mean ft = K and variance cr^ = viK. The response of the detector 
is the sum of signals over all PMTs of the detector. = 10^ events were simulated. 

The MC data were fit with the gaussian response function and with the response 
function based on the generalized gamma- distribution. The results of the fit are pre- 
sented in Tableland FigjGl The mean values and the normalization are reproduced 
well for the gaussian and generalized gamma line shapes; the difference in variances 
is within the statistical precision of the method. The value for the gaussian case 
excludes the hypothesis of the normal line shape; in the case of the non-gaussian 
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shape we have a good match of the data with the model {x^/n.d.f. = 111.6/116, 
the number of degrees of freedom (n.d.f.) here is the number of bins used in the fit 
with the number of free parameters subtracted) . We have found no difference when 
applying method A or B (see Appendix A and B) to estimate of parameters of the 
non-gaussian line shape. 

As it is noted above, Prescott in [6] obtained a precise line shape for the case of 
an exponential single photoelectron response fi{x) = ^e~«, x > 0, it reads: 



where Ji is a modified Bessel function of the first kind for an imaginary argument. 

The slope of an exponential distribution coincides with its mean value, i.e. qi = 
a. The variance of the single electron exponential response doesn't depend on 
parameter a and is vl^^ = 2. It is clear that formula f fTSll can't be directly applied 
to fit the real scintillation shape. The way to solve this problem was pointed out 
in [6]: it is enough to treat a = as a scale parameter, the variance in this case 
will scale as ^/a and the mean value as a. In order to preserve the mean value and 
variance in the original scale, we multiply /i by a scale parameter s = ^ = 
and as before set gi = 1: 



Now formula (fTBl) can be used to fit the scintillation line, the results are pre- 
sented in Table [U Comparing the values one can see that the quality of the 
fit with Prescott formula is worse than in the case of the fit with the general- 
ized gamma function, but much better than in the case of the fit with the normal 
distribution. The quantitative comparison of the models can be performed using 

2 

Fischer's F-distribution gnificance test: ^ = F(a, i/, z/), where z/ is a num- 

Xi 

ber of the degrees of freedom and a is a confidence level [23]. Solving equation 
F{a, 116, 116) = 1883/111.6 with respect to a one can exclude the gaussian shape 
with a c.l. more than 99.99%. The scintillation line shape is described better by 
Prescott's formula (as can be seen from the comparison of values in Table[T]) and 
the exclusion c.l. is smaller, but Prescott's model fails to describe the data with 
high precision as the generalized gaussian distribution does. 

The obtained results have demonstrated very weak sensitivity of the real line 
shape to the shape of the s.e.r., so one can choose any convenient s.e.r. shape in 
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/i 




Norm (xlO^) 


X^/n.d.f. 


MC input 


150.00 


14.18 


1.000 




Gauss 


150.01±0.05 


14.19±0.03 


l.OOOiO.OOl 


1883/116 


Gen. gamma 


150.02±0.05 


14.17±0.03 


l.OOOiO.OOl 


111.6/116 


Prescott 


149.52±0.05 


14.19±0.03 


l.OOOiO.OOl 


329.0/116 



Table 1: Characteristics of three different fits of the monoenergetic line. 



order to invert formula (l5|). 

5.2 ^"^C beta spectrum: MC model of the experimental data 

The major part of the background in the ultra-pure CTF in the energy region up to 
200 keV is induced by /3-activity of ^^C [24], which is present in the organic liquid 
scintillator at the level of 10^^^ g/g. The /3-decay of ^"^C is an allowed ground- 
state to ground-state (0+ —^1^) Gamow- Teller transition with an endpoint energy 
of Eq = 156 keV and half life of 5730 years. The end-point of the decay is used 
in CTF to establish the energy scale, thus the precision of the modeling of ^^C 
spectrum defines the precision of the energy scale calibration. 

The beta energy spectrum with a massless neutrino can be written in the fol- 
lowing form |25| : 

dN{E) ~ F{Z, E)C{E)pE{Q - EfdE (17) 

where 

E and p are the total electron energy and momentum; 

F{E^ Z) is the Fermi function with correction of screening caused by atomic 
electrons; 

C{E) contains departures from the allowed shape. 

For F{E, Z) we have used the function from [26] which agrees with tabulated 
values of the relativistic calculation [27]. A screening correction has been made by 
Rose's method p8] with screening potential Vq = 495 eV. The ^"^C spectrum shape 
factor can be parametrized as C{E) = 1 + aE (see [29] for more details), the value 
of the parameter a was fixed at the value a = —0.7 MeV~^. 

The deviations of the light yield from the linear law have been taken into ac- 
count by using the ionization deficit function /(/cb, E), where is Birks' constant 
[30] . To calculate the ionization quenching effect for the scintillator on the base of 
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Figure 6: Comparison of the MC generated monoenergetic response fit using the 
normal (left) and generalized gamma (right) distributions. To the non-critical eye 
the both fits are comparable in the region /x±2(T, however, the deviations in the tail 
for the gaussian distribution are evident. The — 111-6 value for the generalized 
gamma distribution is close to the number of the degrees of freedom (n.d.f.=116) 
while for the normal distribution the x^/n.d.f — 1883/116 excludes the normal- 
distribution hypothesis. 
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pseudocumene, we used the KB program from the CPC library |3T]. The value of 
the ionization quenching parameter ks = 0.017 cm^^MeV"^ was fixed at the value 
found by independent experiments. The radial dependence of the mean registered 
charge on the point of interaction inside the detector has been accounted for with 
the /^(r) function, obtained from the experimental data (see [3j). For convenience 
the value of the function at the detector's center was assumed to be the unity, 

fR{0) = 1. 

The response of the detector for an event of ^^C decay was generated in the 
following way. First, the event energy E was generated according to the spectrum 
f flTl) . and the position of the event was generated in assumption of uniform dis- 
tribution of ^^C decay events in the detector volume. Then the mean number of 
p.e. has been defined, registered for an event of energy E occurring at distance r 
from the detector center, taking into account detector's non-uniformity and non- 
proportionality of the light yield on the energy: 

Q{E,r) = A-E-fn{r)-f{kB,E), 

where A is the scintillator specific light yield measured in photoelectrons per 
MeV. 

Then in each event for each PMT the mean value of registered number of p.e. 
has been defined, and the registered p.e. number K was generated according to 
the corresponding Poisson distribution. Finally, the registered anode charge was 
simulated by using a gaussian approximation of the PMT signal with mean n = K 
and variance = ViK. The response of the detector is the sum of the signals 
over all PMTs of the detector. = 5 x 10^ event were simulated, that corresponds 
approximately to 3 years of continuous data taking with the CTF detector. 

The exponential underlying background has been added to the ^''C /3-spectrum 
to simulate the realistic situation. We have taken the parameters of the exponential 
observed in the CTF detector. This background is mainly due to the external 7's 
from decays of elements from ^^^U and ^sa^jj^ chains in the water shield. 

5.3 ^^C beta spectrum: fitting MC data with model function 

The real detector response to uniformly distributed events is not spatially uniform. 
To take into account the additional pulse height variance we exploit formula [3]: 
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al = {l+vl)Q + VTQ\ (18) 

where 

Q = A ■ E ■ f{kB, E) ■ fji is the mean total registered charge for the events of the 
energy E uniformly distributed over the detector volume. J'r is the mean 
value of the fnir) function over the detector volume; 

Vi = jf^^ J2i=i ''^ ^i'^u is the relative variance of the PMT single photoelectron 
charge spectrum (fij averaged over all PMTs of the detector {NpM in 
total) taking into account the i-th PMT relative sensitivity Sj. For the 
CTF detector this parameter has been defined with a high precision 
during acceptance tests [8] and turns out to be tJT = 0.34; 

A is the scintillator specific light yield measured in photoelectrons per 

MeV; 

Vt is the relative variance of the photon transfer efficiency, mainly due 

to the spatial non-uniformity of the detector. Among other additional 
contributions there is the intrinsic scintillator line width, the precision 
of the detector calibration, the precision of zero signal definition, etc. 
There is now need to keep these additive parameters apart, so in the 
model we have left the only parameter. In the MC modeling these 
additional contributions were set to zero, but, nevertheless, parameter 
Vt remained free, see discussion below. 

The MC spectrum was modeled with a sum of two components: (1) convolution of 
the ^"^C beta spectrum with the detector resolution function with 3 free parameters: 
total normalization N, light yield A, and additional variance vt; (2) an additional 
exponential background with 2 free parameters. 

The final model function S{Q) has 5 free parameters and is presented as: 

/dE 
N{E{Q'))—Res{Q, Q')dQ' + ExpBkg{Q), (19) 

where Res{Q, Q') is the detector response function, and N(E) is the ^^C beta- 
spectrum (ITTl) . 
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Figure 7: Residual of the fit of the data using the normal and generalized gamma 
distributions (the region up to 100 p.e. is shown). The residual of the fit with the 
normal distribution (upper plot) has two fake peaks in the region of the ^^C tail. 
This is a typical situation for the resolution function mismatch. The fit of the same 
data with the generalized gamma function (lower plot) has no pronounced artifacts 
in the region of the ^^C beta-spectrum tail. 





A 


Norm (xlO^) 


Slope 


X^/n.d.f. 


MC input 


391.8 


5.000 


100.0 




Gauss 


387.8±0.3 (-13(t) 


5.174±0.010 (+17(t) 


99.2±0.5 {-2a) 


279.7/214 


Gen. Gamma 


394.0±0.3 (+7 a) 


5.033±0.008 (+4(7) 


100.0±0.3 (Oa) 


211.3/214 



Table 2: Parameters of the model fitting the CTF MC ^"^C spectrum. Errors cited 
for each parameter are 68% c.l. errors obtained while studying the x^-profile. The 
value in parenthesis near every fitting parameter gives a deviation from the nominal 
value in units of the standard deviation for the corresponding parameter. 
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Figure 8: Fit of MC ^^C spectrum with a model function (only the region up to 100 
p.e. is shown). The fit region 30-250 p.e. corresponds to 91-681 keV. 

The results of the fit of the experimental data with the gaussian and non- 
gaussian line shapes in 30-250 p.e. region, are presented in Table [2l Again, the 

is much better for the non-gaussian line shape. The comparison of the models 
excludes the gaussian shape on the c.l. of 98% (solution of F{a, u, u) = 279.7/211.3 
with u = 214 gives a = 2 x 10"^). 

This time relatively big deviations in parameters have been found when apply- 
ing different resolution functions. The deviations for parameters are bigger than 
statistically allowed, so it should be treated as systematic errors. As it follows from 
Table [21 the error in the light yield definition for the case of the gaussian line shape 
is —1%, the error of the total normalization is +3.5%. With the generalized gamma 
function the error in light yield is smaller: +0.6%, the same error has the total 
normalization. 

It is not implicitly assumed that additional broadening of the scintillation line 
shape {vtQ"^) is distributed in the same way as the main contribution {l+vi)Q. The 
statement is not true in general, especially for big Q values where vtQ'^ term can 
dominate in the response. In our case the main term dominates, that is confirmed 
by the quality of the fit, so the precise distribution for the additional line broadening 
can be neglected. The price paid for this simplification is the observed systematical 
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deviations. 

When fitting the monoenergetic line from a— decays of ^^^Po without selecting 
the detector central region the quality of the fit is much worse at the left side of the 
peak. In the case of ^^C spectrum these imperfections on the left side are covered 
due to the fast decrease in the spectrum and the gaussian shape is justified. On the 
right side the proper description of the scintillation line tail is important because of 
the same fact of the fast decrease of the spectrum. In the case of the monoenergetic 
line the true shape of the distribution of the mean values over the detector volume, 
has to be taken into account. 

6 Conclusions 

An approximation of the real line shape of the scintillation detector with the gen- 
eralized gamma distribution has been proposed. The approximation describes the 
ideal scintillation line shape better than the widely used normal distribution. Two 
parameters of the proposed function are uniquely defined by the first two moments 
of the detector response or by the first two even moments. The computational com- 
plexity of the resolution function calculation is comparable to that of the normal 
resolution. 

It has been demonstrated that the ideal detector response to many photoelec- 
trons {/I > 8) loose the sensitivity to the shape of the single electron response of a 
photomultiplier and the only important parameter is the s.e.r. relative variance. In 
analytical calculations any convenient function can be used instead of a real s.e.r. 

While for the relatively "flat" experimental spectra one can hardly expect the 
enhancement of the overall quality of the fit, in the case of the fast-varying distribu- 
tions, such as tails of the /3— spectrum, the use of the proposed resolution function 
allows one to exclude the artifacts associated with resolution mismatch, and avoid 
systematics errors as demonstrated by the example with the ^^C spectrum fit. 
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Appendix A 



An approximate solution of system (11211 can be obtained using the following expan- 
sion [32j: 

r(a + |) _ / 1 1 5 21 \ 

r(a) ~ V " 8" ^ 128^ + 1024^3 ~ 32768a4 + -J (20) 

For big fi the expansion converges fast because of a ~ /x. Taking three first terms 
and substituting (3 in the first equation, we obtain a simple quadratic equation 

f{a) = 1 - — + 



8a 128«2 + ^2 

with the only positive root: 



2£t 



"° = 16(1-^^) ' ^^'^ 

which gives the solution with a relative precision of ~ lO^'^ for > 10. A more 
accurate solution can be obtained by using more terms from the expansion (l20l) . 
Assuming that more accurate solution has a form a = ao + Aa and developing /(a) 
and two remaining terms from ( l20l l into a Tailor series keeping only a linear term 
with respect to Aa, we obtain a linear equation for Aa with the following solution: 



II — 5aQ 



Aa = J (22) 



Equation (l22|l gives the relative precision of the parameter estimation of < 10 
at /i = 20, at /i = 100 it is ~ 10"^. 



Appendix B 

In radio-engineering the generalized gamma-distribution variant are widely used to 
describe radio waves propagation in fading environment. One of the most popular 
is the m-distribution proposed by Nakagami [16] in the functional form 

1 (m)iZ'" 
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where fl = W^, and m is the inverse of the relative variance of i?^. The advan- 
tages of this equation are simple rules to calculate the parameters. 

In fact, for the even moments of f fTTll the system of two equations for a and j3 
will not contain gamma- functions. Using the parameters a and (3 we can write the 
second and the fourth moments: 



^4 ^ ^-2 r(2+Q) — ^2 . ( r,2 j_ si 

The solution of this system is 



2 « 

(23) 



^ -/3 



(?)' 

a = — 'i 



■ ^^^^ 

In order to use (l24l) . we should require the equivalence of the first two even 
moments of (fTTl) to those of the ideal scintillator response, which can be easily 
calculated with (I9|): 



g4 = /i(l + 6/i + V + ^i(3 + 2/i) + 2v^{3 + 8/i + 2iJ)) + (^q^ 
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